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Analysis of soft optical modes in hexagonal BaTiC^: 
transference of perovskite local distortions 
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We have performed detailed first-principles calculations to determine the eigenvectors of the zone- 
center modes of hexagonal BaTiC>3 and shown that the experimentally relevant low-energy modes 
(including the non-polar instability) can be represented as suitable combinations of basic local polar 
distortions associated with the instability of the cubic perovskite phase. The hexagonal structure 
provides a testing ground for the analysis of the influence of the stacking of TiOg octahedra: the 
occurrence of relatively high-energy chains of dipoles highlights the importance of local effects related 
to the coherent hybridization enhancement between Ti and O ions. Our results provide simple 
heuristic rules which could be useful for the analysis of related compounds. 
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Barium titanate has two structural polymorphs: the 
cubic perovskite type (c-BT) and its hexagonal modifi- 
cation (h-BT) with six formula units per unit cell. While 
c-BT has been one of the best studied ferroelectric ma- 
terials for decades , most of the work on the structural 
and dielectric properties of h-BT is quite recent [0 . As 
shown in Fig. [l| h-BT is also composed of TiOg groups, 
albeit with a different stacking than the perovskite form. 
It seems well established that the hexagonal polymorph 
undergoes two zone-center structural phase transitions: 
at 222 K from the high temperature PQ^/mmc hexag- 
onal phase to a second non-polar C222i phase, and at 
74 K to a ferroelectric P1\ phase. The first transition is 
associated with the softening of an optical mode and the 
second attributed to a shear strain instability; but a de- 
tailed analysis is lacking due to the absence of structural 
information on the two low-symmetry phases, and little 
is known about the microscopic origin of the instabilities. 

On the other hand, the discovery of a giant LO-TO 
splitting in h-BT by Inoue et ai. j| suggested that its 
ferroelectric modes have the same origin as those of c- 
BT. In the cubic phase the ferroelectric instabilities can 
be essentially described as chains of dipoles that originate 
in the movement of Ti ions relative to their surrounding 
06 octahedra with a minor distortion of the later. Addi- 
tional evidence in support of this view is provided by the 
successful use of polar local modes by Zhong et al. in 
the construction of an effective Hamiltonian for c-BT. In 
view of the basic structural similarities between the cu- 
bic and hexagonal forms of BaTi03, it is meaningful to 
ask whether the local modes that describe the unstable 
branches in c-BT can somehow be transferred to h-BT 
and serve as a basis to discuss the low-energy distortions 
of the structure. 

Here we show the results of first-principles calculations 
that provide for the first time structural information on 



the low-symmetry phases of the hexagonal polymorph of 
BaTi03 and reveal the microscopic nature of the modes. 
Our analysis proves that the structure of the experimen- 
tally found zone-center optical soft || modes in h-BT is 
indeed characterized by the same distortions of the TiOg 
octahedra that are relevant in c-BT, leading to similar 
chains of dipoles in the hexagonal structure. 

In h-BT, the experimentally found optical soft modes 
are: the zone-center instability that drives the phase 
transition at 222 K and transforms according to the E2 U 
irreducible representation (irrep) of 6/mmm, and the 
A2 U ferroelectric mode that softens (though remaining 
stable) in the temperature range of the C222i phase and 
is responsible for the giant LO-TO splitting. Our calcu- 
lations agree with this experimental evidence. We per- 
formed a full ab-initio Q relaxation of the thirty-atom 
h-BT structure, resulting in lattice parameters a=10.68 
and c=26.053 a.u. (to be compared with experimental 
values of 10.77 and 26.451 j?J, respectively). The five free 
internal coordinates are also in excellent agreement with 
the experimentally determined ones (within 1%). After 
computing the force-constant matrix at T and diagonal- 
izing it within the subspaces of the appropriate symme- 
tries, we found an unstable Fj2 U mode || and a soft but 
not unstable A 2 „ mode. [It should be noted that if the 
calculations are performed using the experimental lattice 
parameters (i.e., at a larger cell volume), the A2« mode 
is found to be unstable (and the E2„ instability is more 
pronounced). The fact that the A2 U mode is very close 
to being unstable could be particularly relevant for the 
phase transition at 74 K.] 

A first analysis of the eigenvectors for both soft modes 
reveals that the Ba contribution is small (around 10% of 
the total mode norm, compared to 4% in the perovskite 
soft mode). Moreover, we checked that if the Ba ions are 
frozen at their high-symmetry positions the modes are 
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still soft, so we do not consider them the following dis- 
cussion, and focus on the distortions of the TiC>6 groups. 

To make the comparison to the perovskite quantita- 
tive, let us consider the polar deformations of the TiC>6 
groups of c-BT. These are shown in Fig. EL where we 
assume that the Ti ion is located at the origin of co- 
ordinates, so only the displacement patterns of the O 
ions need to be considered. For each spatial direction 
a=x, y, z we have two symmetry-adapted distortions de- 
noted by Si t0l and s 2 ^ a and transforming according to 
the Ti„ (vector like) irrep of m3m, the point group of 
the regular octahedra of c-BT. In terms of this basis, the 
tetragonal ferroelectric distortion in c-BT (along x for 
concreteness) can be written as Q.Q9si x + 0.73s2z, with 
the distorted octahedra exhibiting point group symmetry 
4mm H . In h-BT there are two kinds of octahedra: those 
centered around Ti ions at 2a Wyckoff positions with 3m 
point symmetry (denoted by TiC^l) in Fig. [I]) and those 
arranged around Ti ions at 4/ Wyckoff positions with 3m 
point symmetry (TiOe(2) in the figure) p0| . The octa- 
hedra in the first set are coordinated in the same way as 
in c-BT, i.e., by sharing O ions with six other octahedra. 
Those in the second set are linked to other 1+3 octahedra 
by sharing one O3 face and single O ions respectively. 

Due to the low (as compared to the case of c-BT) sym- 
metry of the octahedra in h-BT, their distortions associ- 
ated to general E 2u and A 2u modes can be decomposed in 
a relatively large number of symmetry-adapted displace- 
ment patterns. Among all the possible ones, we restrict 
ourselves to those of c-BT type and check if they can ac- 
tually account for the structure of the soft modes [ pd| . 
For instance, a general A 2u distortion leads the crys- 
tal to a phase with space group P63mc, in which the 
Ti06(l) groups reduce their point symmetry to 3m (see 
Table ||). As shown in Fig. g.b, c-BT distortions in the 
Six — s iy = s iz component combination (Ti ions move 
towards O3 faces as in the rhombohedric phase of c-BT) 
produce this symmetry breaking. In the case of an E 2u 
distortion, TiOe(l) reduces its point symmetry to 2, and 
the appropriate c-BT mode has the form Si X — —Si V (or- 
thorhombic), as shown in Fig. ||.c. 

For the two soft modes of h-BT, we considered sepa- 
rately the various classes of octahedra, computed from 
our ab-initio eigenvectors the displacement of the O ions 
relative to the Ti ion, and performed a projection of the 
resulting distortion field into the c-BT type symmetry- 
adapted modes. The results (last column of Table |) 
present two main features: First, almost 100% of the to- 
tal structural change associated with both soft modes can 
be described in terms of the c-BT type polar distortions 
(normalization is chosen in such a way that, for instance, 
for the first row in Table § we have (0.62 2 + 0.78 2 ) x 100 = 
99.3%). Second, the components sy, s 2 are always similar 
in magnitude to those of c-BT (0.69; 0.73) and present a 
positive s\/ s 2 ratio, which implies that the Oq octahedral 
cage moves almost rigidly relative to the Ti ion also in 



h-BT. We have also computed the Born effective charge 
associated with the ferroelectric A 2u soft mode and found 
it unusually large (Z* = 11.29), further confirming the 
relation with the (rhombohedric) ferroelectric instability 
of c-BT (for which Z* = 9.956) 

We have proved then that at a local level the soft modes 
in h-BT can be described by the same distortion vec- 
tors that determine the c-BT polar instability. In the 
crystal as a whole, these local polar distortions lead to 
chains of dipolcs, which points at the E 2u instability and 
the softness of the A 2u mode of h-BT as being caused 
by Coulomb destabilizing forces, as it happens in the 
cubic perovskite. The ferroelectric A 2u soft mode, po- 
larized along z' , is roughly depicted in Fig. [l]. In the 
E 2u distortion the chains of dipoles lay on the x'y' plane 
and alternate in orientation with a zero net polarization 
(the resulting C222i phase could be informally consid- 
ered anti- ferroelectric rather than paraelectric). 

From first-principles studies of the c-BT phase it is 
known that parallel dipole chains are very weakly cou- 
pled, so that a transverse modulation of a chain-like in- 
stability is not energetically relevant, and, therefore, un- 
stable TO normal modes exist almost in the whole Bril- 
louin Zone (BZ) |0|| . [The only exception are k points 
near k# = —(1,1,1), for which we have an anti-phase 
modulation of the Ti displacements (Ti=>0<S=Ti-0) in 
the three spatial directions, so the long-range destabi- 
lizing forces are always canceled.] If this view is taken 
to its logical conclusion, we could expect to find more 
zone-center soft modes in h-BT, corresponding to the 
other possible distributions of chains of dipoles. Table || 
enumerates all the possibilities. Apart from the already 
discussed A 2u and E 2u modes, our ab-initio calculations 
show that there is one E\ g mode that is indeed rather 
low in energy, while the ferroelectric E\ u and the E 2g 
modes that are dominated by the movement of Ti ions 
are quite hard. In order to explain this result, let us re- 
mark that for the E 2u and A 2u soft modes the distortion 
is such that if an O ion is approached by one of its two 
Ti neighbors the second Ti ion moves away from it. This 
reflects the hybridization of the Ti 3d and O 2p electronic 
states, which has been shown to be essential for the oc- 
currence of the c-BT ferroelectric instability |l4| . It can 
be checked that any other zone-center arrangement of the 
chains of dipoles results in either two Ti ions approach- 
ing one O ion (for example, if the two Ti ions in one of 
the O3— Ti— O3— Ti— O3 groups depicted in Fig. |l| move 
in the same way in the x'y' plane, there is at least one 
oxygen of the shared face that is approached by both) 
or in the second titanium not moving away from an oxy- 
gen. In the former case (Ei u and E 2g ) the effect of the 
hybridization is lost and the corresponding modes are 
hard. In the latter (E\ g and Bi g ), the hardening is not 
as strong. Thus, we conclude that the particular stack- 
ing of the TiOg groups in h-BT causes (through this local 
effect) the relatively high energy of some chain-like dis- 
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tortions. We can then formulate two "rules of thumb" 
for the characterization of a given locally polar distor- 
tion as low-energy: First (Hi): "There need to be chains 
of dipoles (without regard for their transverse modula- 
tion)". Second (72.2): "The distribution of such chains 
must lead to a coherent hybridization enhancement be- 
tween Ti and O ions, where the word coherent means 
that the destabilizing effect is lost when two Ti ions ap- 
proach the same O ion". These heuristic rules could be 
used to predict the occurrence of locally polar soft modes 
at other k points of the BZ of h-BT, as well as in other 
structures with TiC>6 octahedra as basic building blocks. 

In summary, first-principles calculations of the charac- 
ter of the zone-center modes of hexagonal BaTiC>3 sup- 
port the physically appealing idea that the experimen- 
tally relevant soft modes (including the non-polar insta- 
bility) can be represented as combinations of local polar 
distortions transferred directly from the cubic perovskite 
form of the compound. Our results lead also to heuris- 
tic rules that provide insight into the influence of the 
arrangement of TiC>6 octahedra on the low-energy dy- 
namics of a structure Jl5| ] . 
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TABLE I. Symmetry breakings of the h-BT phase TiO e 
groups associated to the A^v, and E^v. soft modes. The fourth 
column shows the combinations of symmetry-adapted c-BT 
type distortions that are compatible with the symmetry re- 
duction (it applies to both s\ and S2). The last column shows 
the projections of the normalized total distortion of the TiOe 
groups onto the s modes of the second column, in the form 
Si;s2. The i?2u mode removes all symmetry elements from 
the TiOe(2) octahedra, and any combination of s a is possi- 
ble. For this case we have listed in the last column the si 
and S2 projections along each of the three spatial directions 
as well as the modulus. 



FIG. 1. Unit cell of the P6 3 /mmc phase of h-BT. For the 
sake of simplicity, only Ti and O ions are shown. The dis- 
tributions of local dipoles corresponding to the E2 U and A2 U 
soft modes are indicated. 

FIG. 2. Polar (vector like) deformations of the m3m (reg- 
ular) octahedra of c-BT. Only the modes polarized along x 
are indicated. The y and z sets are analogous. 

FIG. 3. Combinations of §j lQ modes that represent the 
symmetry breaking of the octahedra in h-BT by the soft 
modes. Panel a) shows a conveniently oriented octahedron 
which will be assumed to have 3m (resp. 3m) symmetry. 
Oxygen ions are labeled as in Fig. ^ and the cartesian axes 
with origin in the Ti ion are indicated. Panel b) shows the 
s x = s y — s z distortion that breaks the 2-fold axes (resp. 
no symmetry element) and is related to the A% u mode. In 
panel c), a s x = — s y distortion results in a 2 (resp. 1) 
point-symmetry, as appropriate for the Ei-a mode. Only si 
modes are shown, but the same combinations apply to §2- 
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TABLE II. Symbolic description of the possible 
zone-center chains of dipoles in h-BT, classified in terms of 
irreps of 6/mmm. The layers of Ti ions are represented along 
the z direction as in Fig. yj (Ti(l) and Ti(2) refer to Ti ions in 
TiOe(l) and TiOe(2) groups respectively). Arrows indicate 
the orientation of the dipoles (horizontal ones symbolize any 
direction in the x'y' plane), and those set in the same type are 
symmetry related. The superscript / marks the ferroelectric 
modes. Displacement patterns that violate 1Z2 (see text) are 
marked with an asterisk. The bottom line shows the mode 
force constants in atomic units (for B\ g an unambiguous as- 
signment cannot be made). 
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